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[i] The Martian polar night distribution of 1.27 //m (0-0) band emission from O 2 singlet 
delta [C^C'Ag)] is determined from an extensive set of Mars Reconnaissance Orbiter 
(MRO) Compact Reconnaissance Imaging Spectral Mapping (CRISM) limb scans 
observed over a wide range of Mars seasons, high latitudes, local times, and longitudes 
between 2009 and 2011. This polar nightglow reflects meridional transport and winter 
polar descent of atomic oxygen produced from C0 2 photodissociation. A distinct peak in 
1.27 jim nightglow appears prominently over 70-90NS latitudes at 40-60 km altitudes, 
as retrieved for over 100 vertical profiles of 0 2 ( 1 A g ) 1.27 //m volume emission rates 
(VER). We also present the first detection of much (x80 ± 20) weaker 1.58 //m (0-1) band 
emission from Mars 0 2 ( 1 A g ). Co-located polar night CRISM 0 2 (' A g ) and Mars Climate 
Sounder (MCS) (McCleese et al., 2008) temperature profiles are compared to the same 
profiles as simulated by the Laboratoire de Meteorologie Dynamique (LMD) general 
circulation/photochemical model (e.g., Lefevre et ah, 2004). Both standard and interactive 
aerosol LMD simulations (Madeleine et ah, 2011a) underproduce CRISM 0 2 ('A g ) total 
emission rates by 40%, due to inadequate transport of atomic oxygen to the winter polar 
emission regions. Incorporation of interactive cloud radiative forcing on the global 
circulation leads to distinct but insufficient improvements in modeled polar 0 2 ( 1 A g ) and 
temperatures. The observed and modeled anti-correlations between temperatures and 
1.27 jim band VER reflect the temperature dependence of the rate coefficient for 0 2 ( 1 A g ) 
formation, as provided in Roble (1995). 

Citation: Clancy, R. T., et al. (2012), Extensive MRO CRISM observations of 1 .27 im 0 2 airglow in Mars polar night and their 
comparison to MRO MCS temperature profiles and LMD GCM simulations, J. Geophys. Res., 117, E00J10, 
doi: 10. 1029/201 1JE004018. 


1. Introduction 

[ 2 ] The winter polar atmosphere of Mars exhibits dis- 
tinctive compositional, thennal, and dynamical regimes that 
are poorly observed due to the lack of solar illumination, 
which leads to limited infrared and visible measurement 
opportunities. The polar winter (and late fall, early spring) 
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atmosphere is separated from the global atmosphere by a 
polar vortex that limits lower level transport into and out-of 
the polar night region, as dramatically indicated by C0 2 
condensation-related increases in polar night argon [Sprague 
et al, 2007] and CO [Smith et al., 2009] abundances. A 
substantial fraction of this atmospheric C0 2 conversion to 
surface C0 2 ice may occur through the formation of C0 2 
snow clouds [Colaprete and Toon, 2002], In any case, C0 2 
clouds appear prevalent in the polar winter atmosphere 
[Neumann et al., 2003; Hayne et al., 2012], whereas dust 
and water ice aerosols have not been spectrally distin- 
guished from such C0 2 clouds [McCleese et al., 2010]. The 
implication of very pure C0 2 and water ice seasonal cap 
frosts [e.g., Langevin et al. , 2007] further points to extremely 
limited transport of water and dust into the winter polar lower 
atmosphere. At upper levels (above ~60 km) transport into 
the polar winter regions is less restricted, as predicted by 
models [e.g., Haberle et al., 1993; Forget et al., 1999] and 
evidenced by the presence of distinctive temperature 
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Figure 1. A schematic of the vertical and latitudinal distri- 
bution of Mars 0 2 ( 1 A g ) emission (kR/km, color contoured) 
for L s = 95- 1 00° as simulated by the LMD GCM. The merid- 
ional stream function (10 9 kg/s) is superimposed in white 
contour lines. Two primary pathways for 0 2 (' A g ) production 
are the photolysis of ozone over solar illuminated latitudes 
(60S-90N) below 40 km, and the three body recombination 
of atomic oxygen over (mostly) unilluminated fall-winter- 
spring polar latitudes above 40 km altitudes. 

increases associated with the downwelling branch of the 
winter hemisphere Hadley circulation \Conrath et al., 1973; 
Smith et al., 2001; McCleese et al., 2008]. This downwelling 
is also characterized by poleward transport of dayside pho- 
tolysis products and subsequent polar nightglow associated 
with their recombination, first evidenced by Mars Express 
(MEX) observations of nitric oxide (NO) ultraviolet airglow 
in southern polar winter night [Bertaux et al., 2005]. 

[ 3 ] Two recent sets of limb profile measurements from the 
MEX OMEGA and SP1CAM experiments [Bertaux et al., 
2012; Montmessin et al., 2011] and the current MRO 
CRISM investigation further indicate chemical transport 
associated with this same upper atmospheric circulation. 
These measurements reflect the meridional transport of 
atomic oxygen, produced at lower sunlit latitudes, into the 
winter upper atmosphere (altitudes above 80-100 km). 
Entrained in the descending branch of this strong upper level 
Hadley circulation, upper level (enhanced) atomic oxygen 
abundances are transported downward over polar latitudes to 
higher pressure levels where three-body reactions, 

O + O + C0 2 — > 0 2 ( 1 Ag) + C0 2 (1) 

O + 0 2 + C0 2 -» O3 + C0 2 (2) 

produce 0 2 (*A ff ) and 0 3 (Figure 1). These polar night, high 
altitude (40-70 km) layers of 0 2 ( 1 A g ) and 0 3 are very dis- 
tinct in origin and occurrence from Mars lower level 
0 2 (* A g ) and 0 3 abundances, which are associated with local 
photolysis of 0 3 and C0 2 , respectively [e.g., Noxon et al., 
1976; Nair et al, 1994]. Figure 1 schematically distin- 
guishes these two regimes of 0 2 (’Ag-) emission in the Mars 


atmosphere, including their distinct spatial distributions. The 
upper level, winter polar 0 2 ( ! A g ) and 0 3 behaviors are in 
fact more closely related to the production of 0 2 ( ! A g ) in the 
nightside lower thermosphere of Venus [e.g., Crisp et al., 
1996]. The transport related origin of these Mars polar 
night oxygen species provides for a unique window into the 
upper level circulation of the Mars atmosphere. The CRISM 
0 2 (’A s ) observations in particular support extensive spatial 
and seasonal comparisons to general circulation model 
(GCM) simulations of this transport. 

[ 4 ] In the following, we present the extended set of 
CRISM 0 2 (* A g ) measurements obtained in Mars polar night 
at northern and southern high latitudes. These limb profile 
measurements sample a wide range of Mars latitudes, long- 
itudes, and local times; as accumulated from a set of dedi- 
cated CRISM limb observations taken at solar longitudes 
(L s ) of 51°, 64°, 75°, 95°, 134°, 165°, 195°, 265°, and 295° 
over the July 2009 - May 2011 period (Table 1). Inversion of 
the CRISM limb radiances for vertical profiles of volume 
emission rate (VER) yields the detailed spatial and seasonal 
distribution of Mars 0 2 ('A g ) polar nightglow, which is 
compared to model simulations generated by the Laboratoire 
de Meteorologie Dynamique (LMD) GCM with photo- 
chemistry [Lefevre et al., 2004, 2008]. Recent modifications 
to this code [Madeleine et al., 2011a] are considered in the 
context of these model-data comparisons for Mars polar 
0 2 ('Ag) nightglow, as well as contemporaneous/co-located 
temperature profile measurements obtained from the MRO 
Mars Climate Sounder (MCS) [McCleese et al., 2008], In 
particular, the radiative influence of clouds are shown to 
significantly impact winter polar 0 2 ('A g ) and temperature 
distributions. We also report CRISM limb detection of weak 
1.58 /i m band emission associated with Mars polar 0 2 ('A g ) 
nightglow. 

2. CRISM Limb Imaging A^) Observations 

[ 5 ] The MRO CRISM instrument is an imaging spec- 
trometer, designed to obtain moderate spectral (7-15 nm), 
high spatial resolution (15-20 m/pixel) surface spectra for 
Mars over the 0.4-3. 9 pm visible/near-infrared spectral 
region [Murchie et al., 2007]. Although the primary CRISM 
science objectives regard surface compositional variations at 
high spatial resolution, a significant number of atmospheric 


Table 1. CRISM Limb 0 2 (* A g ) Nightglow Observation Summary 


Date 3 

Year b 

MY C 

T d 
L S 

Pole" 

Profiles 1 

' 0 2 (‘A g ), T g 

July 10-11 

2009 

29 

301° 

N 

4 

0 2 (‘A g ) only 

February 10-11 

2010 

30 

50° 

S 

9 

0 2 ( 1 A i; ), T 

April 7 

2010 

30 

74° 

S 

11 

0 2 (‘A ff ), T 

April 28-29 

2010 

30 

84° 

s 

11 

o 2 (‘a s ), t 

May 26 

2010 

30 

96° 

s 

12 

O^'Ajj), T 

August 22-23 

2010 

30 

137° 

s 

23 

0 2 ( Ag). T 

October 17 

2010 

30 

166° 

N.S 

12 

0 2 (‘A g ) only 

December 5-6 

2010 

30 

193° 

N.S 

14 

0 2 ('Ag) only 

March 31 -April 1 

2011 

30 

265° 

N 

4 

0 2 (‘A ? ) only 

May 14-15 

2011 

30 

293° 

N 

4 

0 2 (‘A g ) only 


a UT month, day of CRISM (MCS) observations. 
b UT year of CRISM (MCS) observations. 

c Mars Year of CRISM (MCS) observations, as in Clancy et al. [2000]. 
d Solar longitude of CRISM (MCS) observations. 
e Polar region observed. 

f Number of CRISM limb profile retrievals for 0 2 ( 1 A g ) nightglow. 
8 Co-located MCS temperature profiles compared for SP winter periods. 
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Figure 2. The grand average (see Figure 4) CRISM limb spectrum of polar night 0 2 (' A,,) 1.27 //m limb 
emission (solid histogram line) is compared against the predicted 1.27 pm lineshape (dashed histogram 
line) based upon; the CRISM channel spacing/locations (dotted vertical lines), the CRISM instrumental 
resolution (solid Gaussian line at right, FWHM = 10.7 nm), and the location and relative strength of 
O^'Ag.) emission lines from 2008 HITRAN (vertical lines and box symbols, scale of vertical axis). 
At the top of the figure, CRISM channel numbers and the percentage of integrated 1.27 /mi emission 
contained in each CRISM channel, based on the observed spectrum (data) and the calculated spectrum 
(model), are provided. 


spectral features are incorporated in the CRISM spectral 
range including multiple band absorptions for C0 2 (1.4, 2.0, 
2.7 Atm), H 2 0 (1.35, 1.9, 2.6 /mi), and CO (2.35 /mi). 
Column retrievals for Mars atmospheric water and carbon 
monoxide have been mapped versus latitude and season on 
the basis of these atmospheric band absorptions [ Smith et al., 
2009]. CRISM aerosol analyses have also been conducted 
based upon broad-band absorptions associated with silicate 
and ice particulate compositions [Wolff et al., 2009; M. D. 
Smith et al., Vertical distribution of dust and water ice aero- 
sols from CRISM limb-geometry observations, submitted to 
Journal of Geophysical Research, 2012]. With respect to this 
current CRISM analysis, 0 2 band emission at 1.27 //m is 
present in CRISM spectra, similar in nadir signal character- 
istics to that reported from MEX OMEGA spectral imaging 
observations by Altieri et al. [2009] . That is to say, the coarse 
spectral resolution of CRISM does not reveal the band 
structure of 0 2 ('A g ) emission, but simply captures the inte- 
grated emission. Figure 2 indicates the CRISM channel 
widths against the band structure of the 1.27 pm emission 
from the electronically excited, singlet delta state of molec- 
ular oxygen, O^Ag). 

Oif 1 A g ) ^0 2 (X 3 Z g ,v" = 0) +hw,\= 1 .27 /tm. (3) 


[6] The (unresolved) band structure is indicated by the 
distribution and relative strengths of the individual rotational 
lines, from the 2008 HITRAN data base (vertical lines 


topped by squares, [Rothman et al., 2009]). CRISM spectral 
channels are indicated by dotted vertical lines channel 
numbers are indicated along the top of Figure 2. The instru- 
mental spectral resolution is roughly twice these channel 
widths, as indicated by the Gaussian solid line to the right of 
the band emission (FWHM ~ 10.7 nm). The solid histogram 
spectrum presents the observed spectral profile of 1 .27 /im 
0 2 ( 1 Af emission, based upon an extensive average of 
CRISM limb observations, as discussed subsequently. The 
dashed histogram spectrum presents the modeled spectral 
profile of 1.27 pm 0 2 ( 1 A„) emission, based upon HITRAN 
data and the CRISM spectral parameters. The percentage of 
integrated 1.27 /im band emission sampled within each 
CRISM channel, based upon the observed spectra, is 
provided across the top of Figure 2, labeled Data. The per- 
centage of integrated 1 .27 pm band emission sampled within 
each CRISM channel, based upon the modeled spectral 
weighting, is also provided across the top of Figure 2, labeled 
Model. The minor disagreement between the modeled and 
observed spectral shapes lies well within the uncertainty of 
the CRISM spectral resolution (i.e., 10.7 nm) at this 
wavelength. 

[ 7 ] The primary focus of CRISM 0 2 ( ! Ag) limb observa- 
tions is dayside vertical profiling in support of atmospheric 
photochemistry studies. Indeed, CRISM limb observations 
of 0 2 ( ! A g ) dayglow constitute a unique atmospheric mea- 
surement relevant to Mars ozone photochemistry and 0 2 (* A g ) 
collisional de-excitation rates in a C0 2 atmosphere. Strong 
seasonal and latitudinal variations in the intensity and vertical 
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distribution of Oi(' A„) dayglow are presented in the 2009- 
2011 CRISM limb spectra, indicating distinct orbital and 
seasonal variations associated with atmospheric water vapor 
distribution [ Clancy and Nair, 1996; Lefevre et al., 2004] and 
perhaps heterogeneous chemistry on water ice clouds [Lefevre 
et al., 2008]. Retrievals for dayside profdes of CM’Ag) 
emission will be provided in following analyses in support of 
such photochemical investigations. However, the current 
analysis focuses on CM 1 A,,) emission observed at high alti- 
tudes over fall-winter-spring conditions of polar night. This 
emission presents a new window in the investigation of 
atmospheric circulation. It also presents a less difficult profile 
retrieval problem relative to dayside CK^A^) radiative 
transfer (RT), and so constitutes our initial presentation of 
CRISM profile retrievals for Mars 0 2 ( 1 A„) emission. These 
retrievals for band integrated CM'Ag) polar nightglow from 
CRISM limb imaging observations are developed in the 
following section. Here, we indicate the frequency and spatial 
characteristics of the CRISM limb observations pertinent to 
measurement of Mars CM 1 A^) dayglow and polar nightglow. 

[8] Dedicated CRISM limb observations began in July of 
2009, with the goal of retrieving vertical profiles for Mars 
atmospheric water vapor, CO, CM’Ag), and dust/ice aero- 
sols as a function of Mars season and latitude. The CRISM 
pointing gimbal does not provide sufficient scan motion for 
limb access, such that it is necessary to perform MRO 
spacecraft yaw maneuvers to place the CRISM field-of-view 
(fov) at the atmospheric limb, in the plane of the MRO orbit. 
In addition, the CRISM detector coolers have degraded 
considerably since 2009 to the degree that CRISM limb 
observations constitute a significant requirement on cooler 
operations. Consequently, CRISM limb observations are 
operationally intensive from both mission and experiment 
perspectives. They also, with the exception of MCS obser- 
vations, preclude observations from other MRO investiga- 
tions. In these lights, CRISM limb operations are limited to 
two full orbits roughly every two months (or 30° intervals of 
L s ). This provides pole-to-pole coverage at 6-12° intervals 
of latitude, for two separated Mars longitudes of ~105W 
and 300W (at the equator, corresponding to Tharsis ridge 
and Hellas basin longitudes). Due to spacecraft safing, 
conjunction, and other operational considerations, the actual 
frequency of CRISM limb observations has deviated sig- 
nificantly from the baseline two month cycle. There was also 
an attempt to provide global limb coverage (equivalent to 
one full day of MRO orbits) in August of 23-24, 2010. 
However, the CRISM detector coolers were not able to 
maintain sufficiently cold detector temperatures for a large 
fraction of those orbits. Table 1 indicates the full set of 
2009-2011 limb orbit observations for which CM'Ag) polar 
nightglow is analyzed here. 

[ 9 ] CRISM limb observations are perfonned with the 
spatial dimension of the visible/near-infrared detector arrays 
aligned parallel to the atmospheric horizon. The atmospheric 
limb is scanned in-track with the CRISM gimbal, from limb 
tangents below the surface to altitudes above 120 km. The 
polar CM'Ag) nightglow is observed to be present only over 
40-70 km limb tangent altitudes. The nominal limb resolu- 
tion is of order 50 meters, however the CRISM spectral data 
are lOx binned horizontally (in detector pixels) and verti- 
cally (in image summing) to reduce data rates. This leads to 
vertical intervals of ^500 meters for which 64 spectra are 


returned along the north or south oriented limbs (depending 
on the location of the Sun). We further average the central 40 
of these spectra to obtain a single limb profile of visible-near- 
infrared radiances, in order to obtain maximum signal-to- 
noise ratios. Variations along the ~30 km limb horizon are 
unlikely to yield useful spatial information, given the in-orbit 
~300 km limb tangent paths. 

3. Limb Profile Retrievals for (^('Ag) Volume 
Emission Rates 

[ 10 ] As indicated in Figure 2, CRISM spectra are appro- 
priate to measuring the integrated 0 2 (' A„) band emission at 
1.27 ;im. CRISM nightglow limb analysis is distinct in 
several respects to the dayglow, nadir analyses of Altieri 
et al. [2009] and Fedorova et al. [2006] regarding MEX 
OMEGA and SPICAM 0 2 ('Ag) observations, respectively. 
Obviously, vertical profile retrievals for 0 2 (' A„) emission 
are supported by limb radiance observations. Spectrally 
dependent absorption and scattering contaminations that 
require detailed corrections for nadir dayside spectra may 
be safely neglected in the analysis of these nightglow 
observations. Furthermore, limb path extinction of the 
0 2 (’A g ) nightglow by suspended dust and ice aerosols 
proves unimportant in the aerosol-free (above 20-30 km) 
polar winter atmosphere. Three such Mars polar nightglow 
measurements were obtained from MEX OMEGA limb 
observations, the analyses of which are presented in 
Bertanx et al. [2012]. Here we discuss the spectral, radiative 
transfer, and profile retrieval aspects of inverting a much 
larger set of CRISM limb observations of CM'Ag) polar 
nightglow for vertical profiles of 1.27 /an volume emis- 
sion rates (VER). We also demonstrate the basic characters 
of CRISM limb radiance profiles and spectra for polar 
night CM'Ag) 1.27 /.an band emission, as well as the first 
detection of Mars 0 2 ('Ag) 1.58 /an band emission. 

3.1. Polar Night 0 2 ( 1 A g ) Limb Spectra for 1.27 /xm 
and 1.58 pm Band Emission 

[ 11 ] Figure 3 presents a moderately bright spectrum of 
Mars CM 1 A„) polar nightglow as observed in April of 2010 
( L s = 74°), at a limb tangent altitude of 50 km near 80S 
latitude. The winter polar night presents negligible limb path 
extinction at these altitudes associated with aerosols, as 
indicated by coincident MCS limb profile measurements. 
Radiance units are presented as reflectance (0.001 x I/F) in 
this case to demonstrate the level of the band emission rel- 
ative to surface spectral reflectances CRISM is designed to 
observe ( e.g ., 0.002 versus 0.2). The 0 2 ( 1 A„) band emission 
signal stands well above the spectral noise in the dark-count 
subtracted background of these nightside limb views. 
Detector/filter boundaries appear as spikes near 1.0, 1.65, 
and 2.7 /jm wavelengths and thermal noise increases rapidly 
beyond 2.7 /tm, as described in Murchie et al. [2007]. 
Figure 4 presents a grand average of over 1000 such CRISM 
limb radiance spectra for the 1 . 1 - 1 . 8 /xm spectral window 
(Figure 4, top). This average limb spectrum incorporates all 
limb spectra observed over the 46-55 km limb tangent 
(aeroid) altitude range, for February 2010- May 2011 
polar night limb observations in which 1 .27 /a n polar night- 
glow is detected. In this presentation, the CRISM spectrum 
is provided in modified absolute radiance units (mW/ 
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Wavelength (microns) 


Figure 3. The CRISM visible/near-lR spectrum of polar 
night limb radiance (in I/F x 1000 units) for a limb tangent 
altitude of 50 km (aeroid height) at 79S, 83W on April 7, 
2010 ( L s = 74°). The 1.27 /im emission line of 0 2 ( 1 A g ) 
appears prominently above the background level, which 
presents discontinuities at detector/filter boundaries and 
increased noise/background levels beyond 2.7 /im wave- 
lengths [Murchie et al., 2007]. 


[m 2 -/im-steradians], or 1/1000 of standard CRISM radiance 
values). This grand average indicates that the majority of 
the 1.27 /im CL^Ag.) band emission is contained in four 
CRISM channels centered at 1.27 //in. The expanded scale 
presented in Figure 4 (bottom) indicates a 3er detection of 
extremely weak 1.58 //in band emission from 0 2 (' Ag). the 
first detection of this emission in the Mars atmosphere. 
The predicted ratio for 1.27 /im/1.58 /itn band emission 
(based upon molecular theory and laboratory measurements) 
is somewhat uncertain, but has been observed in the Venus 
nightside lower thermosphere with a ratio of 78 ± 10 
[Piccioni et al., 2009]. The observed Mars band ratio is 
80 ± 20, where the large uncertainty is significantly impacted 
by uncertainty in the spectral baseline. This 0 2 ( 1 Ag) band 
radiance ratio is in agreement with but considerably more 
uncertain than that determined for Venus (where peak 
0 2 (’A g ) emission rates are 5-10 times greater). This result is 
primarily a testament to the dynamic range of CRISM spec- 
tral sensitivity with respect to nightglow detection ( i.e ., in the 
absence of a scattered light continuum level). 

3.2. Polar Night 0 2 ('Ag) Limb Emission Profiles 

[ 12 ] CRISM limb profiles of 0 2 ( 1 A g ) band emission are 
calculated as the difference between a two channel sum 
centered at 1.267 /im (channel numbers 147,148 in Figure 2) 
and the sum of bounding channels centered at 1.257 and 
1.283 /nn (channel numbers 146 and 150). This minimizes 
spectral noise in the retrieved band emission and provides 
optimum subtraction of spectral background radiance asso- 
ciated with issues such as incomplete dark count subtraction. 
Analysis of the CRISM limb average spectrum indicates that 
the central two channel sum incorporates 75% of the inte- 
grated 1.27 /im band emission, versus 10% of the band 


emission contained in the two bounding channels (see 
Figure 2), such that the channel difference effectively 
determines 65% of the total 1.27 /im 0 2 ( ! A g ) band emis- 
sion. By comparison, the predicted proportion of integrated 
1 .27 /.mi band emission for this channel difference is 64%, 
based upon the CRISM spectral resolution (FWHM of 
10.7 nm), channel wavelengths, and 1.27 /im band structure 
as calculated for an atmospheric temperature of 160 K (also 
Figure 2). For band emission calibrations, we adopt the 
observed 65% proportion determined from the average 
CRISM limb spectrum of 1.27 //m band emission (Figures 1 
and 3) and estimate a 2% uncertainty in this calibration fac- 
tor. This uncertainty is also roughly equivalent to the effects 
of a ±20 K perturbation to atmospheric temperature on the 
0 2 ('A ff ) band structure. 

[ 13 ] Such calibrated 1.27 /im emission profiles are con- 
structed for all CRISM polar night limb observations. For 
the purposes of the current analysis, polar night regions are 
defined by latitudes poleward of 70NS during fall-winter- 
spring seasons. Figure 5 presents a set of four such 1.27 /nn 
limb radiance profiles, indicating both peak intensity and 
limb tangent altitude variations. In all presented profiles, the 
altitude scale is calculated with respect to the Mars aeroid 
surface. The limb radiance units for Figure 5 are (Mega) 
Rayleighs, a standard unit of airglow line emission (47t x 
brightness units of 10 6 photons/[cm 2 /second/sterad] [Hunten 
et al., 1956]). These units correspond to previous ground- 
based [e.g., Novak et al., 2002; Kmsnopolsky, 2003] and 
spacecraft-based [ Fedorova et al., 2006; Altieri et al., 2009] 
presentations of Mars 1.27 /im band emission. They also 
support equivalent comparisons to model 1.27 //m volume 


Mars Polar 0 2 ( 1 A ) Nightglow (grand average) 



1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 


Wavelength (/im) 


Mars Polar 0 2 (’A o ) 1.58 /zm Nightglow (grand average) 



Figure 4. (top) A grand average of all CRISM polar night 
limb spectra containing substantial 1.27 /im emission over 
2010-2011 and at limb tangent altitudes between 46 and 
55 km is presented for the 1.1 to 1.9 /im spectral region. 
Standard CRISM radiance units (xlOOO, i.e. mW) are 
employed. The 1.27 /im line emission of 0 2 ('A g ) is well 
characterized in this average, whereas very weak 1.58 /im 
emission is less apparent, (bottom) Expansion of the 1.3— 
1.9 /im spectral region to illuminate the first detection of this 
weak 1.58 /im emission in the Mars atmosphere. The lcr 
absolute uncertainty in Y-axis radiance is 5% (see text). 
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CRISM Limb Radiance Profiles: Polar Night 0 2 ('A g ) 



Figure 5. Four CRISM CL^A^) 1.27 //m limb radiance 
profiles are plotted versus limb tangent aeroid height. These 
limb radiances are presented in units of Mega Rayleighs (see 
text), as calibrated for the total band emission and CRISM 
spectral characteristics. Significant seasonal variation in the 
amplitude and peak altitude of limb emission is indicated 
among the set of high latitude (81-87°) NP ( L s = 165, 
293°) and SP ( L s = 50, 137°) observations. The portion of 
the L s = 165° profile below 40 km exhibits CH^A^) associ- 
ated with low solar elevation photolysis of ozone. 

emission rates (in units of kilo Rayleighs/km). Four periods 
of CRISM limb observations are represented in Figure 5, 
spanning a range of Mars seasons at northern and southern 
high latitudes. Significant spatial (latitude and longitude) 
and local time (LT) variations also occur within each season/ 
observing period, as demonstrated in Figure 6. In this case, 
all CRISM limb observations obtained over 70S-90S lati- 
tudes among CRISM limb observations on April 7, April 28, 
and May 28 ( L s = 74-96°) are distinguished by latitude bins 
of 70S-80S (dashed lines) and 80S-90S (solid lines). 
CRISM polar limb observations extend over the South pole 
(but not the North pole, as yet) to obtain both AM and PM 
local times (LT) of observation. As discussed in following 
model-data comparisons, the apparent larger variability of 
CH^Ag) emission over 70S-80S versus 80S-90S latitude 
bands primarily reflects larger diurnal rather than spatial 
variation at lower latitudes. In addition, all observed periods 
exhibit significant (10’s of %) increases in average (over LT 
and longitude) CLl’Ag) nightglow emission toward higher 
latitudes. 

3.3. Limb to Profile Matrix Inversion 

[ 14 ] The derivation of CL^Ag) volume emission rates 
(VER) from CRISM limb profiles of 1.27 //m radiance 
requires geometric inversion from vertical profiles of inte- 
grated limb path emission to vertical profiles of per unit 
volume emission. Given emission noise levels and the nar- 
row vertical region of C^ 1 Ag) emission, we have adopted a 
non-linear matrix inversion technique constrained to yield 
positive-only solution values [Clancy et al., 1982], The 
radiative transfer (RT) of 1.27 /im limb emission is oblig- 
ingly simple in the Mars polar night regions, due to 


negligible scattering associated with zero or minimal solar 
illumination coupled with the absence of significant molec- 
ular or aerosol opacities at the 40-60 km altitudes of polar 
1.27 /xm nightglow. Upper limits for aerosol extinction in 
the polar night are provided by coincident Mars Climate 
Sounder (MCS) limb profiling, albeit at thermal infrared 
wavelengths. Nevertheless, minimal 1.27 /im aerosol extinc- 
tion along the limb path of 0 2 (' A„) emission is indicated for 
currently measured Mars aerosol particle sizes (R^ of 1-3 /im). 
For conditions of high incident angle solar illumination 
(October and December 2010 observations), visible and near- 
infrared portions of the CRISM spectral limb radiances 
indicate negligible aerosol opacities at 40-60 km altitudes. 
Self-absorption is also negligible (limb path r < 10 4 ) due to 
the very low 0 2 abundances. Consequently, the simplest of 
limb RT conditions are assumed, in which the 0 2 (' A,,) 
radiance is approximated by a vertical profile of 0 2 ( 1 A„) 
1.27 /im VER integrated along a limb path geometry. The 
inversion weighting matrix, P(i, j), is constructed as simple 
limb path weights to the VER profile, E(j). The linear matrix 
equation becomes: 

J2P(i,j)-EU)=S(i), (4) 

j=\,nl 

where S(i) are the limb radiance observations and the sum- 
mation is conducted from layer 1 to the top layer of solution, 
nl. However, P(i, j) values below a given tangent layer, i, of 
observation (i.e., j = 1 to i-1) are zero. In order to specify non- 
negative solution E(j), or VER, we construct a non-linear 
matrix equation in which the solution vector E(j) is expressed 
as e x(j) and the matrix equation becomes: 

Y,P{iJ)-‘? U) =S{i). (5) 

j=l,nl 



CRISM Limb 1.27 /am Radiance (Mega Rayleighs) 

Figure 6. CRISM limb radiance profiles of 1.27 /im 
0 2 (' A„) emission, as in Figure 5, from SP winter limb 
observations obtained in April and May of 2010 ( L s = 74- 
96°). Latitude bands of 70-80S and 80-90S are distinguished 
as dashed and solid lines, respectively. The much greater 
variability apparent at 70-80S latitudes primarily reflects 
LT (diurnal) rather than spatial variation. 
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1.27 fj . m Limb Radiance and Retrieved VER (22Aug 1 0-87S.61 W) 




Figure 7. A non-linear matrix inversion is adopted to invert 
CRISM limb radiance profiles of (left) integrated 0 2 ( 1 A,,) 
1 .27 /tm band emission to (right) vertical profiles of 1 .27 fi m 
volume emission rates (VER in kilo Rayleighs/km). The solid 
line in Figure 7, left, compares the observed 1.27 fim band 
emission for an August 2010 limb observation ( L s = 137°) 
against the RT model limb radiance profile (dashed line) asso- 
ciated with the retrieved 1 .27 /im VER (solid line in Figure 7, 
right). The solution retrieved VER profile is specified for dis- 
crete vertical layers indicated by horizontal lines at the right 
side of Figure 7, right. 


The partial differential matrix, PD(i, j), for X(j) iterative 
solution is simply: 


PD{i,j) — P(i,j)'E(J). (6) 


[ 15 ] We have experimented with a variety of layer thick- 
nesses and numbers. Minimum layer thicknesses of 2 km are 
chosen for the primary altitude region of 0 2 ( 1 A S ) nightglow 
between 46 and 56 km. Increasingly broader layers are 
employed for solution over the lower boundary 30-46 km 
and upper boundary 56-80 km regions. A total of 10 layers 
are determined for the full 30-80 km range of solution. 
Solutions for 12 layers lead to negligible changes in derived 
VER over the altitude region of significant 0 2 ( 1 A g ) 
nightglow. 

[ 16 ] Figure 7 presents the observed (solid line) and fitted 
(dashed line) limb radiance profiles (Figure 7, left), and the 
retrieved VER profile (Figure 7, right) for a CRISM limb 
observation obtained in August 22 of 2010 (L s = 137°). The 
retrieval algorithm adopts an initialized, uniform VER pro- 
file (equivalent to ~15 kRay/km), and quickly converges to 
solution VER profiles that provide accurate fits to all 
observed nightglow limb radiance profiles. The solution 
layers, indicated by horizontal bars on Figure 7 (right, right 
side), provide adequate representation of the VER 1.27 jim 
nightglow profile of VER within the uncertainties of the 
limb radiance observations. Figure 8 presents 0 2 (’A g ) VER 
profile retrievals associated with the four limb radiance 
profiles presented in Figure 5. 

[ 17 ] The uncertainties in retrieved 0 2 (*A ? ) VER for 
these polar night conditions (where aerosol extinction is 


negligible) include contributions from absolute (constant 
scaling for all profiles) and measurement noise (per individ- 
ual retrieved profile) error sources. The absolute la uncer- 
tainty is 5%, as dominated by the radiometric calibration 
error for CRISM radiances in the 1 .27 /im wavelength region 
[. Murchie et al., 2009; D. Humm, personal communication, 
2011]. This includes the small (2%) uncertainty associated 
with calculation of integrated 0 2 ( 1 A s ) band intensity from 
the individual CRISM channel differences. The noise-limited 
uncertainty for individual retrieved 0 2 ('A g ) VER is 4 kR/ 
km, as determined from the binned (vertically and horizon- 
tally) noise level of the CRISM limb radiances and the 
covariance matrix associated with vertical profile retrievals. 

4. C^ 1 A g ) Volume Emission Rate (VER) Profiles 

[is] As a measure of our current understanding of chem- 
istry and dynamics within the polar middle atmosphere, we 
present the full set of retrieved CRISM 0 2 ('A g ) VER pro- 
files in the context of LMD GCM VER profiles and coin- 
cident MCS temperature measurements. The characteristics 
of winter polar temperature and 0 2 ( 1 A g ) nightglow profiles 
are determined in large part by upper level (50-100 km) 
circulation into the Mars polar regions. This circulation is 
globally forced in association with spatially and temporally 
varying distributions of solar heating and thermal cooling. 
Polar winter dynamics is particularly dependent on “distant” 
circulation at low latitudes, as well as upper level boundary 
conditions. As a consequence, the observed behavior of 
polar 0 2 ('Ag) nightglow can only be understood in the 
context of the global Mars circulation. We begin with short 
descriptions of the comparison LMD GCM 0 2 (’A g ) simu- 
lations and MCS temperature profiles. 

4.1. LMD GCM Photochemical Modeling 

[ 19 ] The CRISM 0 2 ( 1 A,,) and MCS temperature profiles 
are compared to those calculated by the LMD GCM with 
interactive photochemistry [ Lefevre et ah, 2004, 2008], 
In the configuration used here, the model is integrated 
on 35 levels extending from the surface up to about 150 km. 



CRISM Limb 1.27 /zm VER (kR/km) 

Figure 8. Vertical profiles of retrieved 0 2 (' A„) 1.27 //in 
VER are presented for the four CRISM limb radiance obser- 
vations presented in Figure 5. 
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and the horizontal resolution is 3.75° latitude x 5.6° longi- 
tude. The photochemical package computes the evolution of 
16 chemical species as well as the O^'Ag.) emission pro- 
duced either by ozone photolysis or by the termolecular 
association of O atoms in reaction (1). For this latter process, 
a net effective yield of 0.75 is adopted for the production of 
0 2 ( 1 A g ). This value proposed by Crisp et al. [1996] is in 
good agreement with the value of 0.7 recently derived for 
Mars and Venus by Krasnopolsky [2010a, 2011]. Both 
studies used a comprehensive scheme of the transfer of 
excitation from the seven electronic states of CL produced in 
reaction (1). Therefore, the production of 0 2 (' A g ) calculated 
by the LMD GCM does take into account the cascade of 
energy and the production of 0 2 ('A^) from the higher 
electronic states of 0 2 . 

[ 20 ] The expression used for the rate coefficient of reac- 
tion (1) is taken from Roble [1995]: 

kl(T) = 9.46 x KL 34 exp(480/r). (7) 

[ 21 ] This exponential law is essentially constrained by the 
laboratory measurements of Campbell and Gray [1973] 
performed at 298 K and 196 K. It predicts a reaction rate 
which is ~20% slower but within the domain of uncertainty 
of the recent measurement of kl carried out at 171 K by 
Smith and Robertson [2008], The strong increase of kl with 
decreasing temperatures is not fully understood from a the- 
oretical point of view, and there are no laboratory data 
available in the domain of T < 170 K explored during the 
Mars polar night. In addition, experimental measurements of 
kl are performed in a bath of N 2 , which is recognized to be a 
less efficient third body than C0 2 . To account for this effect, 
kl and all three-body reactions are multiplied by a factor of 
2.5 in the LMD model. This value has often been adopted in 
the modeling community since the proposition of Nair et al. 
[1994] but has not been precisely detennined for reaction 
(1). In summary, due to the lack of laboratory data relevant 
to the Mars atmosphere, the value of kl used in the LMD 
model simulations (and in all other models) is subject to a 
large uncertainty, especially at low temperatures. 

[ 22 ] Once produced by reaction (1), 0 2 ( 1 A g ) emits a 
photon at 1.27 pm by radiative decay with a characteristic 
lifetime r, or is quenched by collision with C0 2 . In the LMD 
model, t is equal to 4460 s from the measurements of 
Lafferty et al. [1998], whereas the quenching rate is 10 -20 
cm 3 s _1 [Krasnopolsky, 2010b]. At the altitudes considered 
in our study (z > 40 km) the slow quenching rate of 0 2 ( 1 A g ) 
by C0 2 plays a negligible role and is not a source of 
uncertainty. 

[ 23 ] Two versions of the LMD GCM are used for the 
comparison. The first version, referred to as the “ standard ” 
version, is that described in Lefevre et al. [2004]. The so- 
called “ interactive aerosol ” version reflects modifications to 
better represent the radiative impact of atmospheric dust and 
water-ice clouds [Madeleine et al., 2011a], In this new ver- 
sion, the dust layer depth is predicted rather than prescribed, 
as implemented by a “ semi — interactive” dust transport 
scheme [Madeleine, 2011; Madeleine et al., 2011b]. This 
scheme is called semi-interactive because the dust opacity 
profiles are predicted by the model, but scaled so that the 


total column opacity matches the dust opacity observed by 
TES. Note that in this paper, the model is constrained by the 
TES opacity observed during MY26, whereas the CR1SM 
limb 0 2 ('A g ) nightglow observations are acquired during 
MY29 and MY30. The dust particles serve as condensation 
nuclei for water-ice clouds, which are also radiatively active. 
The radiative effect of both dust and water-ice clouds 
depends on the size of the particles, and their radiative 
properties are constantly updated in the GCM as their spatial 
distribution and particle size evolve. 

[ 24 ] In presenting the two distinct model cases, one might 
hope to find confirmation of the specific influences of cloud 
radiative forcing in the CRISM observations of 0 2 ( 1 A g ) 
polar nightglow. However it is also instructive to present 
comparisons for these two models in the context of model 
sensitivity to current uncertainties in Mars global circulation, 
particularly with respect to poleward meridional circulation 
at upper atmospheric levels. 

4.2. Coincident MCS Profiles of Temperature 
(and Aerosols) 

[ 25 ] Winter polar temperature profiles in the 40-60 km 
regions exhibit the same peaked structure as 0 2 ( 1 A g ) night- 
glow. In both cases, this is a direct consequence of polar 
convergent flow from upper levels associated with strong 
Hadley circulation in the fall-winter-spring atmosphere. 
Atmospheric temperatures also directly influence 0 2 ( 1 A,,) 
emission rates through the temperature dependent rate coef- 
ficient for 0 2 (’A ff ) production by three body recombination 
of atomic oxygen (equations (1) and (7)). In these respects, 
measurements of polar winter 0 2 ( 1 A„) and temperature pro- 
files illuminate comparable atmospheric processes. For the 
2010 February-August period of CRISM limb observations, 
we have obtained MCS temperature and aerosol profile 
retrievals that are closest in time and space. While the MRO 
MCS and CRISM instruments do not obtain truly coincident 
limb views, the every-orbit operation of MCS allows fairly 
close coincidence in measurements between these experi- 
ments (typically within two hours of time, 2 degrees of lati- 
tude, and 15 degrees of longitude). The description of MCS 
profile retrievals may be found in Kleinbohl et al. [2009]. 
Key attributes include full global coverages at equatorial 
local times of 3 PM and 3 AM, 5 km vertical resolution over 
0-80 km altitudes, and simultaneous temperature, dust and 
water ice profile solutions. Our use of the co-located MCS 
aerosol profiles is restricted to a determination that dust and 
ice aerosol extinction along the limb path of polar night 
0 2 (’A g ) emission is negligible. With respect to MCS tem- 
perature comparisons, we point out that MCS polar winter 
temperature profiles have already been compared to stan- 
dard LMD GCM temperature profile simulations [Forget 
et al., 1999], These comparisons clearly demonstrate that 
the model does not predict accurate adiabatic heating pro- 
files associated with polar downward convergence 
[McCleese et al., 2008], The combined CRISM 0 2 ('A,,) 
and MCS temperature profile comparisons to LMD GCM 
simulations presented herein provide their first comparison 
to interactive aerosol GCM simulations of polar winter 
temperatures. 
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February 2010, L s =50°: South Polar 0 2 Nightglow 





Figure 9. Averaged Oy'Ag) volume emission profiles (VER) from: (left) all CR1SM retrievals on 
10-1 1 February 2010 ( L s = 50°) for latitude bins of 70-80S (dashed lines) and 80-90S (solid lines); (mid- 
dle) standard LMD GCM simulations co-located and contemporaneous with each CRISM measurement for 
the same latitude bins; (right) interactive aerosol LMD GCM simulations co-located and contemporaneous 
with each CRISM measurement for the same latitude bins in this late southern fall season. 


4.3. CRISM/LMDGCM 0 2 ( 1 A g ) Comparisons 

[ 26 ] We present CRISM profile retrievals of 0 2 ( 1 A g ) 
VER against two separate LMD GCM model simulations, 
employing standard and interactive aerosol models as 
described above. In all cases, the model and comparisons are 
co-located (latitude, longitude) and contemporaneous (same 
L s and LT). We segregate the model-data by Mars season 
into three categories; southern winter, northern winter, and 
northern-southern equinoxes. Given the preponderance of 
2010 observations, the most extensive portion of CRISM limb 
observations pertains to southern winter conditions ( L s = 50- 
137°). All of these southern winter periods correspond to 
CRISM limb observations that extend over the pole to 
obtain the highest latitudes (87S) and substantial diurnal 
coverage (where all hours still remain in polar darkness). 
The August 2010 period of multiple orbit coverage also 
enriches this data set relative to the others. The northern 
winter pole observations were obtained in 2009 ( L s = 301°) 
and 201 1 ( L s = 265, 293°) with limited high latitude (lat < 
82N) and no AM local time coverages. The equinoctial 
periods consist of two 2010 periods ( L s = 166, 193°), 
for which nearly polar night conditions were obtained at 
high latitudes (>77NS) for both poles, with limited diurnal 
coverage. For these model-data comparisons of 0 2 ( 1 A g ) 
VER, we consider the average 0 2 ('Ag) behaviors (vertical 
extent, peak amplitudes); and in the case of the L s = 74-137° 
periods, their global (latitude, longitude) and temporal 
distributions (Ls and LT). We conclude this section with a 


consideration of the key parameters, atmospheric and atomic 
oxygen densities, which contribute to outstanding differ- 
ences between observed and modeled 0 2 (*A g ) polar night 
emission rates. 

4.3.1. The 2010 Southern Winter 

[ 27 ] In Figures 9-11, we compare the CRISM 0 2 ( 1 A„) 
VER (Figures 9-11, left) to simulations for 0 2 ( 1 A g ) VER 
from the standard (Figures 9-11, middle) and interactive 
aerosol (Figures 9-11, right) LMD GCM models; for L s 
periods of 50° (February 10-11, 2010), 74-96° (April 7, 
April 28-29, and May 26, 2010), and 137° (August 22-23, 
2010), respectively. In each case, we show the averages of all 
observed and modeled 0 2 ( 1 A^) VER profiles over the date 
intervals and within two latitude bins of 70S-80S (dashed lines) 
and 80S-90S (solid lines). The late southern spring season of 
February 2010 (L s = 50°, Figure 9) presents the poorest agree- 
ment between observed and modeled 0 2 ( 1 A„) VER profiles 
among all of the observed periods. Both model profiles present 
peak 0 2 ( 1 Ay VER values above 60 km over 80-90S, versus the 
observed peak altitude of 5 1 km. The observations also indicate 
increasing 0 2 ('A g ) VER toward the pole (i.e., larger over 80- 
90S versus 70-80S), which is typical for all of the CRISM polar 
winter observations but most clearly displayed in this period. In 
contrast, the standard model presents fairly constant 0 2 ('A g ) 
VER over 70-90S for this season, and the interactive aerosol 
model presents poleward decreasing 0 2 (' A g ) emission at lower 
altitudes (45-55 km) for this season. The April, May 2010 
period (L s = 74—96° Figure 10) corresponds to southern winter 
solstice and provides perhaps the best agreement between 
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L s =74-96°: South Polar 0 2 Nightglow 
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Figure 10. Averaged 0 2 (*Ag) volume emission profiles (VER) from: (left) all CRISM retrievals on 7, 
28-29 April and 26 May 2010 ( L s = 74-96°) for latitude bins of 70-80S (dashed lines) and 80-90S (solid 
lines); (middle) standard LMD GCM simulations co-located and contemporaneous with each CRISM 
measurement for the same latitude bins; (right) interactive aerosol LMD GCM simulations co-located 
and contemporaneous with each CRISM measurement for the same latitude bins in this southern winter 
solstice season. 


observed and modeled 0 2 ( 1 A„) VER profiles, in terms of alti- 
tude distribution and latitudinal gradient. The interactive aerosol 
model provides slightly improved observational agreement 
for the altitude of peak emission, relative to the standard 
model. Figure 11 presents 0 2 ( 1 A ;; ) VER profiles for late 
southern winter (August 2010, L s = 137°), when both 
models present higher 0 2 ( 1 A„) emission above 55 km than 
present in the CRISM observations. The interactive aerosol 
model provides somewhat better agreement with the obser- 
vations, in terms of both vertical and latitudinal gradients of 
0 2 ('A g ) emission. Overall, these 2010 model-data compar- 
isons for the average 0 2 (*A ff ) emission profile in southern 
polar night indicate modest improvements in agreement for 
interactive aerosol versus standard model simulations. Both 
models tend to bias Oa/Ag) emission to higher altitudes, 
relative to the observations, for late fall and late summer 
periods bounding the southern summer solstice. In contrast, 
the CRISM observations show minimal seasonal variation in 
0 2 (' A„) vertical or latitudinal distributions over this same 
period. Hence, the observations indicate a more constant 
expression of winter polar Hadley circulation over an 
extended seasonal range (L s = 50 = 137°) than exhibited by 
either model. 

[ 28 ] Spatial/diumal variations in the 2010 observations of 
southern polar 0 2 ( 1 A„) nightglow provides somewhat stron- 
ger support to the interactive aerosol versus standard model 


simulations. The average profiles presented in Figures 9- 
1 1 correspond to a range of longitudes, local times (LT), 
and latitudes observed within each latitude bin at each 
observational period. Variations associated with these 
parameters are evident in the limb radiance profiles of 
Figure 6 for the April, May period of 2010 ( L s = 74- 
96°). A significant fraction of this variation corresponds 
to changes in LT of observation as the spacecraft passes 
over the south pole. Such behavior may indicate the 
influence of solar tides through modulation of the vertical 
transport of atomic oxygen at high latitudes. The LT 
coverage of CRISM observations itself is fairly limited and 
somewhat asymmetric about the pole, with irregularly sam- 
pled LT intervals of 2-10 PM and 1-5 AM. Given the sparse 
latitude/LT coverage of the February CRISM observations, 
we focus on the April-August period of southern winter 
( L s = 74-137°). Figures 12 and 13 illustrate the latitude/LT 
variation of O^/Ag) polar nightglow in southern winter for 
the April-May ( L s = 74-95°) and August ( L s = 137°) peri- 
ods, respectively. Both figures present nadir column inte- 
grated values of 0 2 ( 1 A,,) VER from CRISM observations 
(Figures 12 and 13, top) and the LMD standard (Figures 12 and 
13, middle), and interactive aerosol simulations (Figures 12 and 
13, bottom). The plotted 0 2 ( 1 A„) VER correspond to vertical 
columns integrated over the 46-75 km region, and so are 
roughly equivalent to nadir views of polar 0 2 (’Ag) 1.27 jim 
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August 2010, L s =137°: South Polar 0 2 Nightglow 
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Figure 11. Averaged 0 2 ( 1 A„) volume emission profiles (VER) from: (left) all CR1SM retrievals on 
22-23 August 2010 ( L s = 137°) for latitude bins of 70-80S (dashed lines) and 80-90S (solid lines); 
(middle) standard LMD GCM simulations co-located and contemporaneous with each CRISM mea- 
surement for the same latitude bins; and (right) interactive aerosol LMD GCM simulations co-located 
and contemporaneous with each CRISM measurement for the same latitude bins in this late southern 
winter season. 


nightglow. They also provide a comparison of the total polar 
production of 1.27 /./m nightglow associated with meridional 
transport of atomic oxygen, between the observations and the 
models. Over this period, the standard and interactive aerosol 
models produce 40% and 30% less 1.27 /zm emission than the 
observations indicate, respectively. Observed latitudinal, LT 
gradients present in these southern winter polar 0 2 ( 1 A g ) VER 
appear more consistent with the interactive aerosol versus the 
standard LMD GCM gradients for the August period 
(Figure 13) in particular. The latitudinal gradient of 70-90S 
0 2 ( 1 A g ) emission is relatively flat over the dayside (PM) 
hemisphere for the observations (Figure 13, top left) and 
interactive aerosol model (Figure 13, bottom), but distinctly 
decreasing toward the pole for the standard model (Figure 13, 
middle). On the nightside (AM) hemisphere, 0 2 ('A g ) emission 
decreases away from the pole for the observations (Figure 13, 
top) and the interactive aerosol model (Figure 13, bottom), yet 
remains poleward decreasing in the standard model (Figure 13, 
middle). The standard model simulation of south polar 1 .27 gm 
nightglow at L s = 137° (Figure 13) effectively presents a 
minimum at the pole, surrounded by a maximum near 75 S that 
extends over all (measured) local times. In contrast, the obser- 
vations present relatively constant 1 .27 /z m emission up to the 
pole from the dayside (PM) hemisphere, whereupon it 
decreases into the nightside (AM) hemisphere. The interactive 
aerosol (Figure 13, middle) simulated 0 2 ( 1 A„) emission 


behaves similarly, although the PM-to-AM gradient is con- 
siderable smaller than presented in the CRISM observations 
(10-20% versus 30-50%). 

[ 29 ] The distinct latitudinal/LT gradients of polar 1 .27 /im 
nightglow presented by the standard and interactive aerosol 
models suggest distinct polar transport morphologies as 
simulated by these models. The agreement between CRISM 
and interactive aerosol model latitudinal/LT distributions for 
the L s = 137° period in particular would favor the polar 
circulation returned by the interactive aerosol model in this 
season. As demonstrated in a following section, the inter- 
active aerosol model also provides improved agreement with 
respect to MCS temperature profiles at this time. The major 
cloud-radiative feature present in the global Mars atmo- 
sphere at this season is the aphelion cloud belt (ACB). These 
extensive water ice clouds tend to reinforce polar meso- 
spheric convergence and associated adiabatic warming in the 
interactive aerosol model because they enhance the lower 
atmosphere source of thermal tides in the ACB, as found by 
Hinson and Wilson [2004]. Tropical clouds enhance the 
temperature diurnal cycle between the surface and 30 km, 
and thus thermal tides, primarily through their absorption of 
thermal radiation emitted by the surface. This induces strong 
warming during the day and cooling at night. Radiative 
effects of clouds on the thermal structure also affect the 
vertical transport of dust and that further contributes wanning 
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Figure 12. The latitudinal dependences of vertically integrated C^Ag) VER (~ nadir VER, see text) are 
presented for (left) PM and (right) AM local times, for (top) CRISM observations and (middle) standard 
and (bottom) interactive aerosol model simulations for the Apr-May period of CRISM observations (L s = 
74-95°). Dashed lines represent least squares linear fits to the latitudinal/diumal trends. Both observations 
and model C^ 1 Ag) VER are fairly constant across PM latitudes to the pole, but decrease from the pole to 
lower AM latitudes in a fairly monotonic fashion. 


of the northern summer tropical atmosphere between 10 
and 40 km [ Madeleine , 2011; Madeleine et al., 2011b]. In 
terms of dynamics, the enhanced thermal tides tend to drag 
the mean flow toward their phase velocity (—240 m/sec). 
This affects the tropical zonal wind from above the clouds 
into the thermosphere. This, in turn, enhances the meridional 
circulation and thus the convergence of mass and its descent 
above the polar region. Hence, the seasonal, spatial and 
diurnal variations of southern winter C^Ag) polar night- 
glow (and temperatures) are affected by cloud radiative 
forcing, as the distinctions in standard and interactive model 
simulations demonstrate. Of course, we cannot argue that 
cloud radiative forcing is unambiguously identified in the 
CRISM polar night 0 2 ('Ag) observations. Improved obser- 
vational agreement is not consistently demonstrated by the 
interactive aerosol model among the all of the observed 
seasons, and both models seriously underpredict the time 
averaged 0 2 ( 1 A g ) polar nightglow. 

4.3.2. The 2009, 2011 Northern Winter 
[ 30 ] CRISM limb coverage of the northern winter pole is 
much more limited than that obtained for the southern winter 


pole in 2010. It consists of strictly PM local time observa- 
tions obtained in July of 2009 ( L s = 301°) and March/ April 
(L s = 265°), May (L s = 293°) of 201 1. Northern polar winter 
highest latitude (80N-90N) measurements are further limited 
to only two 82N measurements obtained in May of 2011. 
Figure 14 presents the comparison of average northern polar 
winter profiles for latitude bins of 70N-80N (dashed lines) 
and 80N-90N (solid lines) against the standard (Figure 13, 
middle) and interactive aerosol (Figure 13, right) LMD 
GCM simulations for contemporaneous (Ls, LT), co- 
located (latitude, longitude) model averages. Both models 
and the CRISM observations indicate increasing 0 2 ( 1 A,,) 
with increasing latitude within the PM, LT ranges observed. 
However, this behavior is not well determined given the 
limited measurements above 80N at this time. In terms of 
northern versus southern polar winter comparisons 
(Figure 14 versus Figure 10), C^Ag) emission appears at 
lower elevations (relative to the surface aeroid) in the north 
versus the south by ~2-4 km in the observations and models. 
However, the northern relative to southern observed intensi- 
ties of 0 2 ( ! Ag) winter polar emission are observed as fairly 
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Figure 13. The latitudinal dependences of vertically integrated 0 2 ( : A g ) VER (~ nadir VER, see text) are 
presented for (left) PM and (right) AM local times, for (top) CRISM observations and (middle) standard 
and (bottom)interactive aerosol model simulations for the August period of CRISM observations (L s = 
137°). Dashed lines represent least squares linear fits to the latitudinal/diumal trends. Similar to the 
Apr-May (L s = 74-95°) behavior, the observations and interactive aerosol model CM 1 A„) VER are fairly 
constant across PM latitudes to the pole, but decrease from the pole to lower AM latitudes. In contrast, the 
standard model 0 2 ( 1 A g ) VER (Figure 13, middle) exhibit a distinctive, deep minimum in O^'A^) 
emission near the southern winter pole, coupled with diurnal symmetry. 
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comparable (within 10%). With respect to model-data com- 
parisons, the observed northern solstice O^'Ag.) emission is 
significantly stronger than present in either the standard 
model (by 30%) or the interactive aerosol model (by 60%). 
Hence, the models underestimate 0 2 ( 1 A„) polar nightglow 
by a significant margin (^40%) in both northern and south- 
ern polar winters. 

[3 1 ] The distinction in 0 2 ( 1 A g ) integrated emission between 
the two models in this season are unlikely to be associated with 
radiative forcing by low latitude clouds such as the ACB, 
which are limited in extent by the warmer perihelion atmo- 
sphere [e.g., Smith et al., 2001; McCleese et al., 2010]. How- 
ever, the northern polar hood is extensive and exerts 
considerable radiative forcing of the global circulation at this 
time. Furthermore, interactive aerosol simulations overpro- 
duce polar hood (and ACB) optical depths by factors-of-two, 
and so exaggerate its radiative effect in such models 
[ Madeleine et al, 201 la; Haberle et al., 201 1]. Hence, over- 
estimation of poleward upper level circulation modification by 


north polar hood clouds may account in part for the difference 
between the two model 0 2 (' A,,) emission rates. 

4.3.3. The 2010 Northern and Southern Equinox 
[ 32 ] Two sets of late 2010 CRISM limb observations, 
obtained on October 17 ( L s = 166°) and December 5-6 
( L s = 193), provide spring and fall polar conditions at the 
northern and southern poles. These are not strictly polar night 
conditions in that solar illumination is present at low eleva- 
tion angles (~10°) for northern high latitudes (above 75N) in 
October and southern high latitudes (above 75 S) in Decem- 
ber. As indicated in Figures 15 and 16, averaged CRISM 
0 2 (' A„) model VER profiles present lower altitude (below 
45 km) increases that reflect 0 2 (' A„) production by ozone 
photolysis, for northern high latitudes in October and 
southern high latitudes in December. Nevertheless, dis- 
tinctive high altitude 0 2 ( 1 A ff ) emission is present in these 
seasons for both poles, associated with meridional transport 
and recombination of atomic oxygen. Each average 
profile presented in Figures 15 and 16 includes 6 to 7 limb 
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Jul2009,MarMay201 1 , L s =266-301°: North Polar 0 2 Nightglow 





Figure 14. (left) Retrieved O^ 1 A g ) VER profiles are presented as averaged from July 10-11, 2009 (L s = 
301°); March 3 1 -April 1,2011 (L s = 265°); and May 14-15,2011 (L s = 293°) CR1SM limb observations. 
CRISM limb observations at 70-80N (dashed lines) and 80N-90N (solid lines) latitudes are compared to 
LMD GCM O^ 1 A g ) VER model profiles for (middle) standard and (right) interactive aerosol simulations 
averaged for the same L s , LT, latitude, and longitude conditions. 


observations that, for the limited numbers involved, do not 
show clear diurnal variation. Average column 0 2 ( 1 A,,) 
emission rates for these periods more closely approach the 
observed values for the standard (—20%) and interactive 
(+10%) models than for the solsticial periods. 

[ 33 ] Notable symmetry is displayed between the observed 
north (solid lines) and south (dashed lines) polar 0 2 ( 1 A g ) 
profiles (Figures 15 and 16, left) in their variation between 
these L s periods, which are similarly offset before and after 
northern fall/southem spring equinox ( L s = 180°). The peak 
altitude of O^Ag) emission varies by 5-10 km between the 
two polar regions, but in opposite sense between the periods. 
This behavior is not represented in the standard model 
simulation (Figures 15 and 16, middle), yet reasonably well 
represented in the interactive aerosol model. Differences 
between the two model 0 2 (’Ag) emission profiles are most 
striking for north polar latitudes in both periods, and reflect 
large (factors-of-five) differences in simulated atomic oxygen 
densities, as discussed in the following section. Model dif- 
ferences in atmospheric density, as affected by differences 
in polar temperature profiles, are less than 30% for the same 
regions. As for the northern winter period, the primary 
distinction in cloud radiative forcing at this time regards 
waxing northern polar hood clouds of substantial optical 
depth that extend to the pole [Benson et al., 2011]. By 
comparison, waning south polar hood clouds at this time are 


less optically thick and much more limited in latitudinal 
extent [Benson et al., 2010], 

4.3.4. Modeled and Observed Relationships Between 
0 2 ( ' A,,) Atmospheric Density, and Atomic 
Oxygen Density 

[ 34 ] The disagreements between modeled and observed 
polar night O^'Ag) emission are substantial in average 
magnitude (10-60%), vertical dependence, and seasonal 
character. Observed 0 2 (' Ag) polar nightglow emission rates 
exhibit much less seasonal variation in integrated intensity 
and vertical distribution than simulated by either model. 
Here we consider two of the most prominent parameters that 
affect O^Ag) production, atmospheric density and oxygen 
density (equation (1)). Once again, we focus on the southern 
winter period of maximum number CRISM 0 2 ( 1 A„) obser- 
vations as well as coincident MCS density observations. 
Correlations with atmospheric temperature, employing 
coincident MCS temperature profiles, are presented in the 
following section. 0 2 ('A g ) production rates are linearly 
proportional to atmospheric density through the three body 
recombination of atomic oxygen. In Figure 17, we compare 
the relationship between 0 2 ( 1 Ag) and atmospheric density as 
observed (Figure 17, left, CRISM/MCS) and modeled 
(Figure 17, right, for the interactive aerosol model). Both the 
model and the observations present a (noisy) anti-correlation, 
rather than correlation, between atmospheric density and 
0 2 ('Ag) polar nightglow, as averaged over the 50-55 km 
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October 2010, L s =166°: South vs North Polar 0 2 Nightglow 
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Figure 15. (left) Averaged 0 2 ( 1 A g ) VER profiles are presented from October 17, 2010 (L s = 166°) 
CRISM limb observations at high southern (75S-90S, dashed lines) and high northern (75N-90N, solid 
lines), and compared to LMD GCM 0 2 (' A„) VER model profiles for (middle) standard and (right) inter- 
active aerosol simulations averaged for the same L s , LT, latitude, and longitude conditions in this pre 
northern fall equinox, southern spring equinox period. 


altitude region of peak 0 2 (' A ;; ) emission. This somewhat 
surprising result is due to a stronger, offsetting anti- 
correlation between atmospheric and atomic oxygen density in 
the polar winter upper atmosphere, as demonstrated in 
Figure 18. Here, atomic oxygen density is plotted versus 
atmospheric density (Figure 18, left) and 0 2 ('A g ) VER 
(Figure 18, right), as simulated by the interactive aerosol 
model. We note that although 0 2 ( 1 A,,) atmospheric densities, 
and oxygen densities vary between the two model simulations, 
the general relationships presented in Figures 17 and 18 do 
not. Figure 18 (left) indicates that model atomic oxygen den- 
sity is anti-correlated with model atmospheric density. This 
relationship, coupled with the quadratic dependence of 
0 2 ( 1 A ff ) emission with atomic oxygen density, explains the 
modeled (and observed) anti-correlation of polar 0 2 ('A g ) 
nightglow with atmospheric density exhibited in Figure 17, as 
well as the modeled correlation of atomic oxygen density and 
0 2 (’A ff ) nightglow presented in Figure 18 (right). That 
atmospheric density and atomic oxygen density are anti- 
correlated simply reflects the altitude-increasing (density- 
decreasing) profile of polar night atomic oxygen volume 
mixing ratio, as forced by atomic oxygen transport in the 
descending polar branch of the upper level circulation. 

[ 35 ] Figures 17 and 18 effectively identify variable atomic 
oxygen density as the primary driver of spatial and temporal 
variations in 0 2 ( : A g ) polar nightglow. This variable atomic 
oxygen density in turn reflects variable poleward transport 


associated with the upper level meridional circulation. The 
distinctions between standard and interactive aerosol model 
simulations of 0 2 ( ! A g ) polar nightglow reflect the influence 
of cloud radiative forcing on this upper level circulation. 
However, both models seriously underpredict 0 2 ( 1 A s ) 
nightglow emission as observed by CRISM over both winter 
poles, suggesting a more pervasive issue with GCM simu- 
lations of upper atmospheric polar conditions. A stronger, 
vertically deeper, and more extended-in-season poleward 
convergence appears to exist over the upper level polar 
winter atmosphere than presented in current MGCM 
simulations. 

4.4. MCS/LMDGCM Temperature Comparisons 

[36] Comparisons of MCS observed and LMD GCM mod- 
eled temperature profiles for the southern polar winter are 
provided inFigures 1 9 (L s = 74-96'" ) and 20 (L s = 137°). Both 
the model and MCS temperature profiles are selected as con- 
temporaneous ( L s , LT) and co-located (latitude, longitude) 
with CRISM limb observations of 0 2 (’A g ) in April 7, 27-28 
and May 26 (Figure 19), and August 22-23 (Figure 20) of 
2010. While we show that the detailed correlation of 40- 
60 km polar night temperatures and 0 2 (’A g ) emission rates 
does not directly reflect the convergence of meridional cir- 
culation in the polar upper atmosphere (below), the general 
profile behaviors of winter polar temperatures and 0 2 (' A,,) 
emission certainly do. Both exhibit distinctive peaks over 
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December 2010, L s =193°: South vs North Polar 0 2 Nightglow 



0 2 (’A g ) VER (kR/km) 




Figure 16. (left) Averaged 02(' A g ) VER profiles are presented from December 15-16, 2010 (L s = 193°) 
CR1SM limb observations at high southern (75S-90S, dashed lines) and high northern (75N-90N, solid 
lines) latitudes, and compared to LMD GCM 0 2 (' A„) VER model profiles for (middle) standard and 
(right) interactive aerosol simulations averaged for the same L s , LT, latitude, and longitude conditions 
in this post northern fall equinox, southern spring equinox period. 


40-60 km altitudes that correspond to the upper level pole- 
ward transport of atomic oxygen and adiabatic heating asso- 
ciated with the same polar convergent circulation. MGCM 
circulation models have not well reproduced the distribution 
of this polar heating, based upon initial LMD GCM com- 
parisons with MCS polar winter temperature profile mea- 
surements [. McCleese et al., 2008]. The presented LMD 
CGM and CRISM southern polar winter 0 2 (' A ;; ) distribu- 
tions also indicate significant model-data discrepancies in 
standard model simulations that are somewhat mitigated by 
incorporation of active-aerosol processes in the LMD GCM 
simulations (notably Figures 13-15). Figures 19 and 20 show 
that, to varying degrees, incorporation of active-aerosol 
processes in the LMD GCM simulations also leads to general 
improvement in comparisons between MCS and LMD GCM 
southern winter polar temperature profiles. This improve- 
ment is not particularly notable in Figure 19 for model/data 
comparisons over April-May, 2010 (L s = 74-96°), in which 
the interactive aerosol (dashed lines in Figure 19, bottom) 
simulated temperature profiles show modestly improved 
agreement with MCS profiles (solid lines) relative to standard 
model profiles (dashed lines in Figure 19, top) for the 70S- 
80S latitude range (Figure 19, left) in particular. More sig- 
nificant improvements in model/data comparisons are 
obtained for the August 2010 ( L s = 137°) period presented 
in Figure 20. In this case, the interactive aerosol simulated 
temperature profiles (dashed lines in Figure 20, bottom) 
compare significantly better with observed MCS profiles 
(solid lines) than the standard model simulated profiles 


(dashed lines in Figure 20, top) for both 70S-80S 
(Figure 20, left) and 80S-90S (Figure 20, right) latitude 
ranges. Hence, the incorporation of interactive aerosol pro- 
cesses in the LMD GCM improves model simulations of 

0 2 (’A g ) vs Atm Density: SP Winter (L g =74-137°, Lot = 70-90S, Alt=50-55km) 




Oj('a g ) VER (KR/km) 0 2 ('A g ) VER (KR/km) 

Figure 17. (left) Observed atmospheric densities (from 
MCS) versus 0 2 ( 1 A,,) emission rates (from CRISM), aver- 
aged over the 50-55 km altitude region from April, May, 
and August ( L s = 74-137°) periods of observation, (right) 
The same interactive aerosol model quantities, selected as 
coincident with observed locations (latitude, longitude) and 
times ( L s , LT). Similar anti-correlations between atmo- 
spheric density and 0 2 ( 1 A g ) emission rates are presented 
by the observations and model. This counterintuitive behav- 
ior is associated with anti-correlation between atmospheric 
density and oxygen density (Figure 18, see text). 
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LMDGCM (interactive), SP Winter (L s =74-137°, Lat=70-90S, Alt = 50-55km) 
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Figure 18. Oxygen densities from the interactive aerosol 
model, selected to be coincident with April, May, and 
August ( L s = 74-137°) CRISM observations, are averaged 
over the 50-55 km altitude region and plotted versus (left) 
averaged model atmospheric densities and (right) model 
0 2 ( ' A„) emission rates. The anti-correlation of atmospheric 
density with oxygen density reflects the strong altitude- 
increasing gradient of atomic oxygen in the polar night 
upper atmosphere. The roughly linear dependence of 
0 2 ( 1 A g ) emission with atomic oxygen density (Figure 18, 
right) reflects the quadratic dependence of O^'A^) produc- 
tion with atomic oxygen density (equation (1)) countered 
by the anti-correlation of atmospheric density (also in 
equation (1)) with atomic oxygen density (Figure 18, left). 


southern polar winter temperatures and 0 2 ( 1 A g ) nightglow, 
although both models still depart significantly from the 
CRISM 0 2 (' A„) and MCS temperature profile 
measurements. 

4.4.1. The Temperature Dependence of Winter Polar 
Emission 

[ 37 ] A characteristic of 0 2 ( 1 A g .) nightglow in the Venus 
lower thermosphere is the correlation of 1 .27 //m emission 
intensity with 10-20 K temperature increases over several 
thousand km regions, and the extreme temporal and spatial 
variabilities of this emission [ Crisp et al . , 1 996; Bailey etal., 
2008; Ohtsuki et al., 2008; Hueso et al., 2008]. Vertical 
downwelling rates associated with such Venus nightside 
regions of intense 1.27 /un emission (VER~500 kR/km) and 
20 K compressional temperature increases are estimated of 
order 20 cm/s [ Bailey et al., 2008]. LMD GCM vertical 
downwelling velocities within the Mars polar winter region 
of 1.27 ftm emission are also of order 10's cm/s. However, 
the pressure level of the Mars polar night peak limb emis- 
sion (50 km, or ^0.006 mbar) is typically ten times smaller 
than that of the Venus nightside lower thermospheric 
emission (^0.06 mbar, assuming a 96 km altitude from 
Drossart et al. [2007], as are peak 0 2 (*A„) emission rates. 
Furthermore, the Mars polar winter 0 2 ( r A g ) emission is 
characterized by similarly reduced temporal and spatial 
variabilities. This particular difference apparently leads to 
very different correlations between local temperatures and 


MCS (solid) vs standard LMDGCM (dashed) for AprMay 2010 (L s =74-96°) 
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MCS (solid) vs interactive LMDGCM (dashed) for AprMay 2010 (L s =74-96°) 
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Figure 19. MCS observed (solid lines) and LMDGCM model (dashed lines) southern winter ( L s = 74- 
96°) polar temperature profiles are compared for (top) standard and (bottom) interactive aerosol simula- 
tions (see text). Model and observed values are binned for (left) 70S-80S and (right) 80S-90S latitude 
ranges. The model and observed profiles are contemporaneous ( L s , LT), co-located (latitude, longitude) 
with CRISM limb observations of 0 2 ( 1 A„) obtained on April 7, 28-29 and May 26 of 2010. The interac- 
tive aerosol model profiles (Figure 19, bottom) provide moderately improved agreement with the MCS 
observed profiles over 70S-80S (Figure 19, left). 
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MCS (solid) vs standard LMDGCM 
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(dashed) for August 2010 (L s =137°) 
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MCS (solid) vs interactive LMDGCM (dashed) for August 2010 (L s =137°) 
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Figure 20. MCS observed (solid lines) and LMDGCM model (dashed lines) southern winter (L s = 137°) 
polar temperature profiles are compared for (top) standard and (bottom) interactive aerosol simulations 
(see text). Model and observed values are binned for (left) 70S-80S and (right) 80S-90S latitude ranges. 
The model and observed profiles are contemporaneous (L s , LT), co-located (latitude, longitude) with 
CRISM limb observations of 0 2 ( 1 A„) obtained on August 22-23 of 2010. The interactive aerosol model 
profiles (Figure 20, bottom) provide significantly improved agreement with the MCS observed profiles 
over 70S-80S (Figure 20, left) and 80S-90S (Figure 20, right) latitude ranges. 


1 .27 /jm VER in the Mars winter polar middle atmosphere 
versus the Venus nightside lower thermosphere. 

[38] In Figure 2 1 , we present plots of Mars southern polar 
winter (74-87S, L s = 74-137°) 1.27 /im VER versus atmo- 
spheric temperature, from CRISM and MCS measurements 
(Figure 21, left) and from the LMD GCM simulation 
(Figure 21, right). Instead of the strong correlation in tem- 
peratures and 0 2 (* Ag) emission observed for Venus, we see 
a noisy anti-correlation for Mars polar night temperatures 
and Oo( 1 Ag) emission. This behavior is presented in both the 
CRISM/MCS observations and the LMD GCM model, to a 
very consistent degree (Figure 21 employs the interactive 
aerosol LMD calculations, but similar results apply from the 
standard model). The dashed lines plotted in Figure 21 
represent least squares fits for exp(/f/7) dependence 
between 1.27 /.an VER and temperatures, with essentially 
identical derived coefficients, A, for the observations (380 ± 
20) and model (430 ± 40). The exp(4/7) functional form 
reflects the temperature dependence of the rate coefficient 
for 0 2 ( 1 A ? ) formation from atomic oxygen (equation (1)). 
The LMD GCM uses the temperature dependence employed 
by Roble [1995]; exp(480/7)) for the Mars 0 2 ( 1 A g ) simula- 
tions (see prior discussion in the LMD GCM description). 
The agreement between the observations and model dem- 
onstrated in Figure 9 provides some support for this 
assumption. The fitted observation and model values for this 
coefficient are perhaps slightly smaller ( i.e ., 380-430 vs 
480), but it does not appear that there is significant offsetting 


positive correlation between 0 2 ( 1 A„) VER and tempera- 
tures. Hence, the correlation of such adiabatic heating with 
0 2 ('A g ) production is not clearly identified and so must be 
weak compared to conditions in the Venus nightside 
thermosphere. 

[ 39 ] As a test of the model temperature-0 2 ('A^) sensitiv- 
ity, we have conducted LMD GCM simulations (interactive 
aerosol) employing the Smith and Robertson [2008] rate 
coefficient for reaction 1 in which the temperature sensi- 
tivity scales as exp(720/7). Employing the same model 
temperature-O^ 1 A g ) analysis above, we derive a value of 
700 ± 50 for A. This is equivalent to the temperature sen- 
sitivity of the Smith and Robertson [2008] rate constant, 
further suggesting this rate constant dominates the model 
tcmpcrature-0 2 ( 1 A K ) correlations. It also indicates that the 
temperature dependence of the Smith and Robertson [2008] 
rate constant leads to disagreement with the CRISM/MCS 
observations, whereas the temperature dependence of the 
adopted rate coefficient from Roble [1995] provides agreement. 

5. Conclusions 

[ 40 ] Vertically (46-75 km) integrated, latitudinally (70- 
90NS) averaged 0 2 ( 1 A„) emission rates of 500-600 kR are 
retrieved from CRISM winter polar nightglow observations. 
The magnitude of LMD simulated polar nightglow is gen- 
erally ~40% lower than these observed rates, for both 
northern and southern winter polar conditions. Simulated 
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0 2 ( 1 A g ) vs Temperoture: SP Winter (L s =74-137°, Lot=70-90S, Alt=50-60km) 
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Figure 21. The correspondence between SP winter tem- 
peratures and Oaf 1 A^) 1.27 jam VER, as averaged over the 
50-60 km altitude (aeroid) region for the CR1SM locations 
(latitude = 70-90S, all longitudes) and times ( L s , LT) of limb 
observations from April-to-August, 2010 (L s = 74-137°). 
(left) An observed anti-correlation between co-located 
CRISM 0 2 ( ' A s ) and MCS temperature observations, (right) 
Very similar anti-correlation between LMD GCM 02( 1 A g ) 
and temperature simulations. The dashed lines in both panels 
indicate a functional fit to the presented 0 2 (' A^/tempera- 
ture distributions, assuming the exp(A/T) dependence of 
the rate coefficient for winter polar (^('A,,) formation 
(equation (1)). Very similar fits are found for the observa- 
tions and model, where the derived coefficients are 380- 
430 K, as compared to the adopted model temperature coef- 
ficient of 480 K or a more recent laboratory detennination of 
720 K by Smith and Robertson [2008] (see text). 
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vertical profiles of winter polar 0 2 (' A g ) emission also show 
much greater seasonal ( L s = 50-137, 266-301°) and vertical 
(46-75 km) variation than presented by observed 0 2 ( 1 A g ) 
profiles. Observed winter 02( 1 A g ) emission profiles exhibit 
more narrowly peaked emission regions centered at 50 ± 
3 km for all of these periods, whereas simulated 0 2 ( 1 A g ) 
emission extends more broadly to altitudes above 60 km 
for the L s = 50, 137° periods bounding winter solstice. 
Vertically integrated 0 2 cA g ) emission rates for two equi- 
noctial periods ( L s = 166, 193°) exhibit better (10-20%) 
agreement between the observed and simulated polar 
regions (700-800 kR). Both equinoctial periods present 
distinctly offset 0 2 ( 1 A g ) emission peaks in the northern 
versus southern polar regions, behavior which is more clearly 
reproduced in the interactive aerosol versus standard LMD 
simulations. The LMD interactive aerosol model also better 
represents the latitudinal and diurnal 0 2 ( 1 A g ) variations pre- 
sented by CRISM measurements as well as MCS polar tem- 
perature measurements at L s = 137°. This southern winter 
period corresponds to the presence of the aphelion cloud 
belt, which leads to strong modulation on Mars thermal tides 
[e.g., Hinson and Wilson, 2004]. Such tidal influence further 
leads to enhanced zonal winds from the upper clouds levels 
(~40 km) into the thermosphere, for interactive aerosol LMD 
simulations [ Madeleine et al., 2011a, 2011b]. This in turn 
increases and diumally modulates meridional circulation and 


mass convergence at southern winter latitudes, and thus mass 
convergence/descent over the pole. Northern polar hood 
clouds may play a related role in the distinctions between the 
equinoctial (^(’A^) profiles simulated by the standard and 
interactive aerosol models. 

[ 41 ] Model simulations and CRISM/MCS observations 
indicate strong anti-correlation between (^(’Ag) emission 
and atmospheric density at 50-55 km level of peak C^ 1 A g ) 
emission, counterintuitive to the relationship expected from 
the production reaction (1) for 0 2 (' A„). This behavior 
instead reflects the sharply increasing mixing ratio of atomic 
oxygen with altitude (i.e., decreasing with density) over this 
region, associated with meridional transport and upper level 
convergence of atomic oxygen in the winter polar atmo- 
sphere. Both model simulations exhibit much stronger cor- 
relations between atomic oxygen density and 0 2 ( 1 A K ) 
emission. Hence, the transport efficiency of atomic oxygen 
plays the most important role in characterizing inadequate 
simulated versus observed winter polar 0 2 (’Ag) emission 
rates. We also show that modeled and observed anti- 
correlations in 0 2 ( ' A„) emission rate and atmospheric tem- 
perature are in close agreement, and mimic the temperature 
dependence of the atomic oxygen recombination coefficient, 
specifically that given by Roble [1995]. This agreement does 
not hold when the Smith and Robertson [2008] rate coefficient 
temperature dependence is employed. In either case, there is 
no indication of correlated temperature and 0 2 ( ! A g ) variations 
as exhibited in the Venus lower thermosphere [e.g., Bailey 
et al., 2008]. 

[ 42 ] We have concluded that both the standard and the 
interactive LMD GCM models fail to provide sufficient 
transport of atomic oxygen to the winter polar regions of 
observed 0 2 (’Ag) nightglow. Perhaps part of this failure is 
related to the outstanding disagreements in observed and 
simulated cloud optical depths associated with the interac- 
tive aerosol simulations [. Madeleine et al. , 2011a; Haberle 
et al., 2011]. It must be noted, however, that polar night 
upper atmospheric circulation, characterized here by the 
adiabatic warming and 0 2 (' Ag) nightglow, is among the 
most sensitive phenomena in general circulation modeling. 
This sensitivity is particularly striking in model inter- 
comparisons [e.g., Mischna and Wilson, 2008]. For instance, 
Forget et al. [1999] compared two models sharing the same 
physical parameterizations and only differing by the numerical 
methods to solve the atmospheric dynamic equations. The 
simulations were usually quite similar, except in the fall and 
winter polar night upper atmosphere where differences up to 
30 K were found [see Forget et al., 1999, Figure 8]. Similarly, 
any changes in poorly constrained processes, such as gravity 
wave drag, horizontal dissipation, convection and turbulence 
(which act near the surface, but influence the tides) have been 
found to significantly affect this part of the atmosphere, often 
without imparting significant changes elsewhere. The 
assumed upper atmosphere conditions can also play a role. 
For instance, in the version of the LMD GCM including a full 
thermosphere, the solar cycle variations of the extreme UV 
flux, which heats the thermosphere above 120 km, can change 
atmospheric temperatures by several Kelvins above 60 km in 
the polar night [Gonzalez-Galindo et al., 2009]. Polar night 
mesosphere observations of tracers such as 02 ('A g ) and 
temperature thus provide unique constraints to test our 
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understanding of the Martian circulation and of planetary 
atmosphere dynamics in general. 

[ 43 ] Last, extensive averaging of polar night CRISM limb 
spectra (over altitude, latitude, and time) reveal the 1.58 /im 
spectral emission of CLf 1 A K ) in the v = 1 vibrational state. 
The observed ratio of C^A^) 1.27 /im to 1.58 /im emission 
is 80 ± 20, consistent with the same ratio observed in the 
Venus nightside lower thermosphere [Piccioni et al., 2009]. 
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